Eight new MSPs in NGC 6440 and NGC 6441 
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Abstract. We report the discovery of five new millisecond pulsars in the globular cluster NGC 6440 and three new ones in 
NGC 6441; each cluster has one previously known pulsar. Four of the new pulsars are found in binary systems. One of the 
new pulsars, PSR J1748-2021B in NGC 6440, is notable for its eccentric {e = 0.57) and wide (Pi, = 20.5 days) orbit. If 
the rate of advance of periastron is due solely to general relativity, we can estimate of the total mass of this binary system: 
2.92±0.20M,;i. This would imply an anomalously large mass for this pulsar, which could introduce important constraints in 
the study of the equation of state for cold neutron matter. 
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DISCOVERY 

The recent discovery of 30 millisecond pulsars (MSPs) 
in Terzan 5 Ol 01, made using the GBT's S-band re- 
ceiver and the pulsar spigot |3] shows that this observ- 
ing system (henceforth GBT/S/PS) has an unprecedented 
sensitivity to MSPs outside the Arecibo sky. This is 
only partly due to the factor of 3 higher gain of the 
GBT as compared to the Parkes telescope. The higher 
observing frequency (~2GHz), large observing band- 
width (600 MHz), and relatively fine frequency resolu- 
tion (~0.78MHz) provide enough frequency and time 
resolution to detect faint MSPs at high dispersion mea- 
sures (DM). An impressive demonstration of this is 
the discovery of PSR J1748-2446ad, the fastest spin- 
ning pulsar known [P = 1.396ms). This was found 
in Terzan 5, where the average DM of the pulsars is 
239cm^^pc. This motivated us to search for MSPs in 
other promising GCs with this system. 

For the first phase of the survey, we selected the globu- 
lar clusters (GCs) NGC 6388, NGC 6440 and NGC 6441 
based on their high stellar interaction rates at the core, 

Despite the large r^ s of these GCs, only two pulsars 
were known in them: PSR B1745-20A (NGC6440A, 
lH) and PSR J1750-37A (NGC 6441A, il). The spin 
periods and DMs of NGC 6440 A and NGC 6441 A are 
288.6 ms and 1 11. 6 ms and 220 and 234 cm^^ pc, respec- 
tively. We assumed that the lack of other known pulsars 
in these high-DM clusters was largely due to a bias of 
the earlier surveys against the detection of MSPs. The 
GBT/S/PS system can greatly reduce this bias. 



With full track observations made in 2005, we discov- 
ered eight new MSPs, five of them in NGC 6440, three 
in NGC 6441 and none in NGC 6388. Preliminary re- 
sults on these objects were presented in Begin (2006)|'7'|. 
Their pulse profiles, together with those of the previously 
known objects, are presented in Figure [T] These pulse 
profiles were obtained by adding all detections of the 
pulsars, carefully excluding all sub-integrations that are 
corrupted by RFI. 

LUMINOSITY DISTRIBUTION 

In Figure |2] we make cumulative plots of the pseudo lu- 
minosities of the pulsars in NGC 6440, NGC 6441 and 
Terzan 5, the latter based on the flux density estimates re- 
ported in 1 1 ] . This allows a direct comparison of the sizes 
of the pulsar populations of these clusters. The observ- 
ing system and observing frequency are the same and the 
DMs of the pulsars are very similar, making the search 
similarly sensitive to fast-spinning pulsars. The software 
used to do the searches is the same, and the software used 
to make the flux density measurements is similar (but not 
exactly the same). All these factors make the flux den- 
sities directly comparable; the pseudo-luminosities will 
then be directly comparable as well if the cluster dis- 
tances are accurate. We have used the latest distance esti- 
mates for these GCs (5.5 kpc for Terzan 5,jl8|], 8.2 kpc for 
NGC 6440 and 13.5 kpc for NGC 6441 0\) to estimate 
the pulsar pseudo-luminosities, but caution that these dis- 
tances still have large systematic uncertainties. 
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FIGURE 1. Average 1950 MHz pulse profiles for the 10 pulsars known in the GCs NGC 6440 and NGC 6441. These profiles 
cover one full rotation. The pulsars in binary systems have an asterisk after their names, and the newly discovered pulsars have 
their names in boldface. The horizontal width of the rectangles indicates the system's total time resolution, including the effects of 
dispersive smearing, relative to each pulsar's spin period. The dip in power after the main pulse of NGC 6440A might be a SPIGOT 
artifact. 
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FIGURE 2. Cumulative distribution of pseudo luminosities 
at 1950 MHz for the pulsars in NGC 6440, NGC 6441 and 
Terzan 5, assuming the most recent distance estimates. 



We have fitted a luminosity function of the type 
to the observed pulsar pseudo- 
luminosity distribution of each cluster, where A^(> L) is 
the number of pulsars above a pseudo luminosity of L 
and ki is the number of pulsars brighter than 1 mJy kpc^. 
We obtain ki = 10.6,7.5 and 7.3 for Terzan 5, NGC 6440 
and NGC 6441 respectively. The clusters seem to have 
similarly sized pulsar populations, as expected from 
their stellar interaction rates: Terzan 5, NGC 6440 and 
NGC 6441, r,. have 6.3, 6.3 and 8.6% of the total stel- 
lar interaction rate of all Galactic GCs. We detect more 
pulsars in Terzan 5 because that cluster is closer to us. 



A SUPER-MASSIVE NEUTRON STAR? 

One of the most interesting new discoveries is 
PSR J1748-2021B. This is a 16.7-ms pulsar in a 
wide {Pj, = 20.55 days) and eccentric (e = 0.57) orbit 
with a ^ O.lMp, companion. This orbital eccentricity 
has allowed a measurement of the rate of advance of 
periastron: d) = 0.00391(18)°yr"'. 

Assuming that co is fully relativistic (as argued in [4]), 
that implies a total system mass of 2.92 ± 0.20Mq. In 
Figure [3] we present the mass constraints on the compo- 
nents graphically, with probability distribution functions 
for the mass of the pulsar and the companion presented 
above and on the right respectively 

The total mass of this system suggests that it is a dou- 
ble neutron star (DNS) binary. However, the probability 
that it lies within the range of neutron star masses that 
have been precisely measured to date, indicated in Fig.|3] 
by grey bars, is only 0.10%. This eventuality would re- 
quire improbably low orbital inclinations, between ~ 4° 
and - 5°. 

The median pulsar mass is 2.74 Mq, with the lower 
and upper 1 a limits at 2.52 and 2.95 Mq. The pul- 
sar has only a 1.0% probability of being less massive 
than 2.01 Mq, and a 1.0% probability of being more 
massive than 3.24 Mq,. If this high mass value is con- 
firmed, it would be the largest pulsar mass measured, 
surpassing that for PSR B1516H-02B in the GC M5 
(1.96+U ?^M0, (T^). Confirmation of such a large neu- 
tron star mass would have profound consequences for the 
study of dense matter, since it excludes almost all known 
equations of state [TJ]. Since DNS systems have simi- 
lar masses, this might even imply that the end products 
of the coalescence of DNS systems might themselves be 
stable as super-massive neutron stars! 
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FIGURE 3. Constraints on the masses of NGC 6440B and its companion. The hatched region is excluded by knowledge of the 
mass function and by sin / < 1. The diagonal dashed lines correspond to a total system mass that causes a general-relativistic d) 
equal or within l-CT of the measured value. The five solid curves indicate constant orbital inclinations. The gray bars indicate the 
range of precisely measured neutron star masses, this goes from --^1.20 Mc-, for the companion of PSR J1756— 2251| 1 1] to 1.44 Mq 
for PSR B 1913+16 |12]. We also display the probability density function for the mass of the pulsar (top) and the mass of the 
companion (right), and mark the respective medians with vertical (horizontal) lines. 
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